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Abstract. We present a detailed phenomenological study of photoproduction of two hadrons, both with 
high transverse momentum, in longitudinally polarized lepton-nucleon collisions. We consistently include 
"direct" and "resolved" photon contributions and examine the sensitivity of the relevant spin asymmetries 
to the gluon polarization in the nucleon and to the completely unknown parton content of circularly 
polarized photons. Our results are relevant for the COMPASS and Hermes fixed-target experiments as well 
as for a possible future polarized lepton-proton collider like eRHIC at BNL. So far, all studies are limited 
to the lowest order approximation of QCD. 
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1 Motivation and Introduction 

The fundamental question of how the spin of the proton is 
composed of the spins and orbital angular momenta of its 
constituents, quarks and gluons, still remains unanswered. 
Over the past 25 years, a series of polarized deep-inelastic 
scattering (DIS) experiments has revealed that the quark 
spins contribute remarkably little to the nucleon spin 1.. 
Measuring Ag{x,^), the spin-dependent gluon distribu- 
tion in the nucleon, in an as large as possible range of mo- 
mentum fractions x, is the prime goal of all current exper- 
iments with polarized beams and targets. In the light-cone 
gauge the first moment of Ag{x, /i), i.e. ^5(2^, ^i)dx, can 
be interpreted as the gluon spin contribution to the nu- 
cleon spin at a momentum scale /i The missing piece, 
the orbital angular momenta of quarks and gluons, might 
be accessible in exclusive processes, but precise measure- 
ments are challenging and rather distant at this point. 

The advent of the Relativistic Heavy Ion Collider 
(RHIC) at Brookhaven National Laboratory (BNL), has 
opened up unequaled possibilities to access Ag over a 
broad range in a; in a variety of high-transverse momen- 
tum, "high-pT" , processes such as, for example, inclusive 
hadron or jet, prompt photon, and heavy flavor produc- 
tion P). In each case, the gluon density prominently con- 
tributes through gluon-gluon fusion and quark-giuon scat- 
tering processes already at the lowest order (LO) approxi- 
mation of QCD. Center-of-mass system (c.m.s.) energies of 
up to ^fS = 500 GeV guarantee that the standard frame- 
work of pcrturbative QCD (pQCD) can be used reliably to 
learn about all aspects of helicity-dependent parton den- 
sities at RHIC. A series of unpolarized "benchmark" mea- 
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surements at RHIC has nicely confirmed the applicability 
of pQCD methods 0] and are the foundation for similar, 
ongoing measurements with polarization. First, very re- 
cent results from the Phenix and Star collaborations at 
RHIC 5 indicate that large and positive gluon distribu- 
tions are disfavored in the range of momentum fractions x, 
0.03 X < 0.2, dominantly probed in these experiments. 
Future, more precise measurements will extend the range 
in X and further close in on Ag. 

The gluon polarization can be accessed also in low en- 
ergy fixed-target experiments like COMPASS Pl at CERN 
and Hermes iZj at DESY. Here one scatters a beam of 
longitudinally polarized leptons off loiigitudinally polar- 
ized nucleons at c.m.s. energies of ~ 18 GeV and 
y/S ~ 7.5 GeV, respectively. Compared to RHIC, Ag{x, fx) 
is probed in a more limited x-range, 0.1 ^ x < 0.2, but 
at smaller momentum scales fi which makes results com- 
plementary. In particular, high-pr hadron pairs, both in 
photoproduction and in deep-inelastic electroproduction, 
have been identified to be the most promising processes 
for a determination of Ag at the low energies available at 
fixed-target experiments [5]. First results for double-spin 
asymmetries are available from Hermes 1^1 (for all pho- 
ton virtualities) , Smc TU^ (electroproduction), and, most 
recently. Compass jllj (photoproduction). These data are 
consistent with only moderate gluon polarizations as we 
will discuss in detail later. 

However, at the smaller c.m.s. energies of fixed-target 
experiments it is much less obvious that standard pQCD 
methods are applicable as straightforwardly as at collider 
energies to analyze data. In fact, experimental results for 
high-pT processes in, e.g., hadron-hadron fixed-target scat- 
tering have been a serious challenge in the past for the 
standard factorized framework where the pcrturbative se- 
ries is truncated at a given fixed order in the strong cou- 
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pling and where possible power corrections are assumed 
to be negligible ■ It is therefore crucial to demonstrate 
first that standard pQCD methods can be used to learn 
about the parton and/or spin content of nucleons in a 
given kinematical regime, for instance, by analyzing un- 
polarized cross sections. Otherwise conclusions about, e.g., 
the gluon polarization might be incorrect. 

In this paper we present a detailed phenomenological 
study of photoproduction of hadron pairs at LO accuracy 
of QCD. Quasi- real photons have the advantage of yielding 
much higher rates than deeply-inelastic electroproduction 
of hadrons. The price to pay is the more involved the- 
oretical framework for photoproduction, where so-called 
"direct" and "resolved" photons contribute to the cross 
section as depicted in Figs. 1 (a) and (b), respectively. In 
(a) the photon simply interacts as an elementary pointlike 
particle, whereas in the latter case the photon "resolves" 
into its parton content prior to the hard QCD interaction, 
for instance, by fluctuating into a vector meson with the 
same quantum numbers. Here, cross section estimates re- 
quire knowledge of the parton content of circularly polar- 
ized photons which is lacking completely at the moment. 
We will demonstrate below, that this does not seriously 
limit the usefulness of this process for the kinematical re- 
gion specific to both Compass and Hermes. Since the 
theoretical framework for hadron-pair production is signif- 
icantly more complex than for single-inclusive cross sec- 
tions, next-to-leading order (NLO) QCD corrections are 
still not complete in the polarized case at the moment "TT. . 
We note that there had been earlier studies of hadron-pair 
production |14| in the wake of first experimental results 
from Hermes However, the phenomenological results 
presented in 14 cannot be compared easily with the re- 
cent (and upcoming) results from COMPASS JT] we are 
aiming at. 

Investigations of interactions between polarized lep- 
tons and nucleons will, hopefully, continue to play a vi- 
tal role in spin physics also in the future. A polarized 
lepton-nucleon collider such as the eRHIC project at BNL 
jl5| . which is currently under discussion, would be the 
next logical step. Besides the unique possibility to access 
Ag{x,^) down to a; ~ 10^'^ in studies of scaling viola- 
tions in DIS, photoproduction processes are particularly 



interesting also at collider energies ^Hl- As we shall show 
below, they will allow to probe for the first time differ- 
ent models for the parton content of circularly polarized 
photons. 

The paper is organized as follows: in Sec. 2 we briefly 
recall the theoretical framework for photoproduction of 
hadron pairs. Section 3 is devoted to detailed phenomeno- 
logical studies. Here we mainly focus on the COMPASS and 
Hermes experiments, where first data are already avail- 
able. We present results for spin asymmetries and discuss 
their possible sensitivity to Ag but also focus on predic- 
tions for unpolarized reference or "benchmark" cross sec- 
tions which allow to probe the validity of the pQCD frame- 
work at low c.m.s. energies and transverse momenta px- 
We include all relevant experimental cuts in our calcula- 
tions. We close this section by studying the prospects for 
hadron-pair photoproduction at eRHIC. We briefly con- 
clude in Sec. 4. 

2 Technical Framework 

We consider the spin-dependent inclusive photoproduc- 
tion cross section for the process 

lipi) N{pn) ^ I'M H,{p,) Hdipd) X , (1) 

where a longitudinally polarized Icpton beam I scatters off 
a longitudinally polarized nucleon target N producing two 
observed hadrons He and Hd in the final state. The pi de- 
note the four-momenta of the particles. Both hadrons He 
and Hd are assumed to have high transverse momenta pt,c 
and pT,d, respectively, ensuring large momentum trans- 
fer. Invoking the factorization theorem |17j we may then 
write the differential cross section as a convolution of non- 
perturbative parton distribution and fragmentation func- 
tions and partonic hard scattering cross sections, 

dAa = ^[d(7++ - rf(T+_] = (2) 

J dXadXb dzc dZd Af{Xa,Hf) ^f'^{xb,tJ'f) 

abed 

X dAa"''^'''^^ {S, Xa,Xi„Pc/ ZcPd/ Zd, fJ.f, f^'f, fJ-r) 

X D^^{z,,^^'f)D^-{zd,^i'f) . (3) 
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In ((21 the subscripts "++" and "H — " denote the helici- 
ties of the colliding leptons and nucleons. S is the total 
c.m.s. energy squared available, i.e., S = (pi + pnY ■ The 
sum in Eq. Q runs over all possible partonic channels 
ah cd with dAa"''^^'^'^ the associated spin-dependent 
LO partonic hard scattering cross sections. The latter can 
be calculated in pQCD order-by-order in the strong cou- 
pling as{^r), with /ir denoting the renormalization scale. 

The Af^{xb, /i/) are the usual spin-dependent parton 
distributions of the nucleon 



(4) 



evolved to a factorization scale /i/, with Xb the momen- 
tum fraction of the nucleon carried by the parton /. The 
subscript + [— ] in Eq. Q indicates that the parton's spin 
is aligned [anti-aligned] to the spin of the parent nucleon. 
The other non-perturbative functions 13^ ^ '^ (zc,d, M/) de- 
scribe the collinear fragmentation of the partons c and d 
into the observed hadrons He and H^, respectively, with 
Zc,d the fraction of the parton's momentum carried by the 
produced hadron. fi'f denotes the final-state factorization 
scale which can be different from fif. 

The experimentally measured cross section for is 
the sum of the so-called "direct" and "resolved" photon 
contributions, cf. Figs. ^ (a) and (b), respectively. 



dAa = dAadir + dAa-^, 



(5) 



We shall note that neither dAaair nor dAa^cs are measur- 
able individually. Both, dAa^ir and dAa^es, can be cast 
into the form of Eq. ^ by defining the parton distri- 
bution functions for a lepton, Af\xa, Hf), appropriately. 
Most generally, they can be written as convolutions. 



with 



I^^AP,i{y)Ar (^x. 



■r-—^^'f 



(6) 



2tt 



2mfyi 



1 i-v. 



(7) 



being the spin-dependent Weizsacker- Williams equivalent 
photon spectrum [T^ that describes the collinear emission 
of a quasi-real photon with momentum fraction y and vir- 
tuality less than some (small) upper limit Qmax off a lep- 
ton of mass mi. Qmax is determined by the experimental 
conditions. 

The explicit form of the polarized photon structure 
function Af'^{x'y,fi) in Eq.® depends on the specifics 
of the interaction that the quasi-real photon undergoes 
in the hard scattering with the nucleon. In the "direct" 
case, depicted in Fig. 1 (a), parton a in ||2Jl has to be 
identified with an elementary photon and hence Xa with 
the momentum fraction y of the photon w.r.t. the parent 
lepton, i.e.. 



in Eq. . If the photon resolves into its hadronic structure 
before the hard scattering takes place, the Af^ in Eq. © 
represent the parton densities of a circularly polarized 
photon. The latter are defined in complete analogy to the 
ones for a nucleon target in Eq. Q). Unlike hadronic par- 
ton distributions, photonic densities consist of a pertur- 
batively calculable "pointlike" contribution, which dom- 
inates their behavior at large momentum fractions x^, 
and a non-perturbative "hadronic" contribution dominat- 
ing in the low-to-mid x-y region. Nothing is known about 
the latter, such that we have to invoke some model for it 
in our calculations below. This will become important in 
the discussion of the numerical results in the remainder 
of the paper. We will demonstrate that measurements at 
low c.m.s. energies, i.e., at Compass and Hermes, are 
to a large extent not affected by the actual details of the 
model. At higher c.m.s. energies, like at a future polarized 
ep collider, one of the physics goals would be a first deter- 
mination of the partonic structure of circularly polarized 
photons. 

Finally, the experimentally relevant double-spin asym- 
metry All is defined as 



dAa — (i(T+_ 



da 



da4 



da 4 



(9) 



AP{x^,^i)^5{l-x^) 



(8) 



The required spin-averaged cross section da in Eq. (|5J) is 
straightforwardly obtained from Eqs. ©-Q by replacing 
all polarized quantities by their appropriate unpolarized 
counterparts. 



3 Phenomenological Applications 

In our phenomenological studies based on the framework 
laid out in Eqs. ©-0 we concentrate on the production 
of pairs of charged hadrons made of light quark flavors. 
In fact, we sum over pions, kaons, and (anti-)protons and 
use the fragmentation functions of KKP 19 throughout 
to model hadronization. All our results will be differential 
in the transverse momentum pT,c of hadron He and inte- 
grated over all kinematically and experimentally allowed 
transverse momenta PT,d of hadron Hd and pseudo-rapi- 
dities ric,d unless stated otherwise. The pseudo-rapidities 
of the hadrons are measured w.r.t. the direction of the 
incident lepton beam. 

For the unpolarized parton densities of the nucleon 
and photon we adopt the LO CTEQ6L [H] and GRV EU 
sets, respectively. To study the sensitivity to the unknown 
gluon polarization of the nucleon we use four different sets 
of spin-dependent parton distributions emerging from the 
GRSV analysis Qfl^- These sets span a rather large range 
of gluon densities Ag all very much consistent with present 
DIS data. Apart from our default "standard" set of GRSV 
with a moderately large, positive Ag, the three other sets 
"Z\(7 = g input", "Z\(7 = input", and "Z\(7 = —g input" 
are characterized by a large positive, a vanishing, and a 
large negative gluon polarization, respectively, at the in- 
put scale of the evolution. 
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Fig. 2. Scale dependence of the unpolarized (a) and polarized (b) photoproduction cross section for two charged hadrons at LO for 
Compass with ^max = 70mrad. The scales in Eq. Q are varied simultaneously in the range {Pt.c +PT,d)/4 < 4:{Pt,c + Pt d) 



as indicated by the shaded bands. The solid lines correspond to the default choice where fi 
impact of lifting the experimental cut p^^c + Pt d > 2.5 GeV^. 
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The unknown parton densities of circularly polarized 
photons arc estimated with the help of two extreme mod- 
els ,23, based on maximal, Af'^{x, /io) — P{x, ^q), or min- 
imal, Af^{x,iJ,o) = 0, saturation of the positivity bound 
at the starting scale /j,o for the evolution to scales /i > /io- 
Both models result in very different parton distributions 
at small-to-medium Xj while they almost coincide 
as x-y —f I due to the dominance of the perturbatively 
calculable "pointhke" contribution in this region. The use 
of the "maximal" set will be implicitly understood except 
when we study the sensitivity of the photoproduction cross 
sections and spin asymmetries to the details of the non- 
perturbative hadronic input to the evolution of Af^. 

Unless stated otherwise, all factorization and renor- 
malization scales, /i/, fi'j:, and fir, in Eq- ® are set equal 
to ^2 ^ /i2 = ^2 ^ = p-^^ + 

3.1 Two-Hadron Production at COMPASS 

With the present setup, the Compass experiment scatters 
polarized muons with a beam energy of E"^ = 160 GeV off 
the deuterons in a polarized ^LiD solid-state target corre- 
sponding to a c.m.s. energy of ~ 18 GeV. On average 
the beam polarization is ~ 76%, and about — 50% 
of the deuterons can be polarized with an average polar- 
ization of -Pd ~ 50% 0. 

Hadrons can be detected if their scattering angle is less 
than ^niax — 70 mrad in the laboratory frame. This accep- 
tance was recently upgraded to ^max = 180 mrad for all 
future runs. In the event selection for a "high-py" sample, 
the charged hadrons have to pass further cuts m] : the in- 
variant mass m{Hc, Hd) of the two produced hadrons has 
to be larger than 1.5 GeV and the sum of the transverse 
momenta squared must exceed Pt c ~^ Pr d > 2. 5 GeV^ 
with both pt^c and px.d larger than 0.7 GeV. In addition. 



the fractions Zc,d of the parent parton's momenta carried 
by the detected hadrons Hc^d are chosen to be z^^d > 0.1. 

The maximal virtuality of the quasi-real photons in 
Eq. lO is taken to be Q^^^ = 0.5 GeV^. The fraction y of 
the lepton's momentum taken by the photon is restricted 
to be in the range 0.1 < y < 0.9. We note that the often 
omitted non-logarithmic pieces in Eq. Q result in a small 
but non-negligible contribution for muons. 

Figure[21shows the dependence of both the unpolarized 
and polarized LO photoproduction cross section Q on 
the unphysical factorization/renormalization scales varied 

in the range (pf.^^ + p'r.d)/'^ < < ^Pt.c + PT,d)- Both 
cross sections exhibit a very large scale dependence which 
is, however, not uncommon for LO estimates. Sets of po- 
larized parton densities with a moderate gluon polariza- 
tion like GRSV "standard" result in an almost vanishing 
cross section as the two "direct" subprocesses, photon- 
gluon-fusion and QCD-Compton scattering, cancel each 
other almost entirely, see also Fig. ^ and the discussions 
below. Even the sign of the polarized cross section can- 
not be determined within the scale uncertainty here. As 
always, the computation of the relevant NLO QCD cor- 
rections is a mandatory task as theoretical uncertainties 
associated with the residual scale dependence are due to 
the truncation of the perturbative series at a given order 
and are expected to decrease significantly beyond the LO 
approximation. Such a calculation for hadron-pair produc- 
tion is a formidable task and still not complete at present 

m 

From a recent calculation of the NLO QCD correc- 
tions to spin-dependent, single-inclusive photoproduction 
of high-pr hadrons we know though that theoreti- 
cal scale uncertainties are only marginally reduced, if at 
all, at the NLO level at c.m.s. energies relevant for COM- 
PASS and Hermes. This is in sharp contrast to what one 
generally expects and what indeed happens at collider en- 
ergies, see, e.g., |25II26| . In addition, the NLO corrections 
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Fig. 3. Double-spin asymmetry All, Eq. at LO for differ- 
ent gluon polarizations in the nucleon (see text). The upper panel 
shows the spin asymmetry with (solid lines) and without (dashed 
lines) including the "resolved" contribution. In the lower panel 
the dependence of All on the two extreme photon scenarios, 
"maximal" (solid lines) and "minimal" saturation (dashed lines), 
is shown. The "error bars" indicate the estimated statistical accu- 
racy for such a measurement at Compass in certain bins of pt,c 
based on an integrated luminosity of Ifb^^. 



in [23 turned out to be sizable and different for the polar- 
ized and unpolarized cross sections thus affecting also the 
spin asymmetries defined in Eq. Q. Clearly, this under- 
pins the delicacy of perturbative calculations in the fixed- 
target regime, i.e., for small c.m.s. energies and transverse 
momenta of only a few GeV. Before drawing any con- 
clusions about the hadronic spin structure, in particular 
Ag, from measurements of spin asymmetries All one has 
to demonstrate the applicability of pQCD methods first. 
For this purpose, an important "benchmark" would be 
the comparison of the relevant unpolarized cross sections 
with theoretical expectations, e.g., given in Fig.|21(a). Un- 
fortunately, such kind of information is still lacking from 
Compass for the time being. Let us stress, that a possible 
discrepancy between experiment and theory at moderate 
c.m.s. energies would not necessarily imply that standard 
pQCD methods are beyond remedy. It would only call for 
further improvements by resumming large terms in the 
perturbative series, for instance, threshold logarithms, to 
all orders in the strong coupling. This is known to lead 
often to a much improved agreement between data and 
pQCD calculations, see, e.g., fI7\ . 



The dotted lines in Fig. El correspond to the unpolar- 
ized (a) and polarized (b) cross sections computed without 
imposing the experimental cut on ^+p^ ^ > 2.5 GeV . 
This cut, which ensures hard scattering, is responsible for 
the cusp observed at around pt,c — 1-4 GeV and for the 
significant reduction of the cross section for smaller pt,c- 

Figure 01 shows our expectations for the double-spin 
asymmetry All, Eq. @, at LO based on the cross sec- 
tions shown in Fig.[21for the default choice of scales. Apart 
from the "standard" set of GRSV polarized parton den- 
sities 22 , we also use the three other sets, introduced at 
the beginning of Sec. 3, with very different assumptions 
about the gluon polarization. In the upper panel of Fig.O 
we study the importance of the "resolved" photon contri- 
bution to the photoproduction cross section By com- 
paring the experimentally relevant spin asymmetry for the 
sum of "direct" and "resolved" contributions (solid lines) 
with All computed for the "direct" part alone (dashed 
lines) one can infer that irrespective of the chosen Ag 
the "resolved" part is non-negligible. It leads to a signifi- 
cant shift in the absolute value of the spin asymmetry and 
neglecting it in the analysis would clearly lead to wrong 
conclusions about Ag. 

The impact of the unknown, non-perturbative parton 
content of the circularly polarized photon on All is ex- 
amined in the lower panel of Fig. (31 by making use of the 
two extreme models j23| also introduced at the beginning 
of Sec. 3. As can be seen, the actual choice of the model 
barely affects the results for the spin asymmetry. Any dif- 
ference between the two results diminishes further towards 
larger pt_c- This can be readily understood as the photonic 
parton densities are probed on average at medium-to-large 
momentum fractions Xj. In this region the partonic con- 
tent of the photon is dominated by the "pointlike" con- 
tribution which is independent of the details of the un- 
known non-perturbative input (2^1 ■ Certainly, this finding 
somewhat simplifies the theoretical analysis of the spin 
asymmetry in terms oi Ag. 

Figure|31also demonstrates that the double-spin asym- 
metry All is sensitive to different model assumptions for 
the gluon polarization. A large and positive gluon po- 
larization yields a sizable negative asymmetry whereas a 
large and negative Ag leads to a positive asymmetry. For 
the "standard" set of GRSV or when Ag = is imposed at 
the input scale of the evolution we find asymmetries close 
to zero. To judge whether a measurement of All can be 
actually turned into a constraint on Ag we estimate the 
expected statistical accuracy 6All for Compass in cer- 
tain bins of pt,c, calculated from 



Here, fTf,i„ denotes the unpolarized cross section integrated 
over the bin considered and £ the integrated luminosity 
for which we assume £ = 1 fb~^. All other quantities are 
as specified at the beginning of Sec. 3. Clearly, the region 
of 1 < pt,c < 2 GeV is the most promising one to obtain 
information about the gluon polarization. At higher pt,c's 



6 



C. Hendlmeier et al.: Hadron-Pair Photoproduction in Longitudinally Polarized Lepton-Nucleon Collisions 



. . . . 1 . . 


11 1 . Ill, 11 11 1 

dir tot 










: dAa^es/dAa^ot 


1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 _ 

Yg : 










Y 

qg - 






T - 

gq . 






T 

qq ■ 

Y 

gg - 








. . . . 1 . . 


, . 1 . . . . 1 . . , , 1 


yq : 



0.5 



1.5 2 

Pt,c [GeV] 



2.5 



Fig. 4. The upper panel shows the "direct" and "resolved" con- 
tribution to the experimentally relevant polarized cross section. 
The lower panel shows the fractional contributions of the differ- 
ent partonic LO channels a + b^c+d to the photoproduction 
cross section. 



the achievable statistical precision deteriorates as the cross 
section drops steeply with pt,c- 

First experimental results [111 find that the spin asym- 
metry (integrated also over pT,c for the time being) is close 
to zero. From this measurement a value of Ag/g = 0.024± 
0.089(stat.) ± 0.057(sys.) at cc = 0.095 + 0.08 (-0.04) and 
scale /i^ = 3 GeV^ was extracted with the help of "puri- 
ties", i.e., a "signal-to-background" separation based on 
Monte-Carlo simulations (llj. Compared to our theoreti- 
cal expectations in Fig.Olthe COMPASS result of All — 
is consistent with a moderate gluon polarization like 
in the GRSV "standard" set. However, as discussed above, 
one should take this result with a grain of salt unless an 
unpolarized reference cross section becomes available from 
Compass to verify the applicability of pQCD methods. In 
particular, one has to exclude that the observed smallness 
of the spin asymmetry is due to the presence of large non- 
perturbative effects. If these are spin-independent they 
would naturally lead to All — irrespective of Ag. 

Next we turn to a closer analysis of how the results in 
Figures [51 and 121 can be understood. To this end we study 
the different contributions to the polarized photoproduc- 
tion cross section separately as illustrated in Figure^l The 
upper and lower panel show the fractional contributions 
of dZ\crdir a-nd dAa^^s and of the different partonic LO 
channels a + b —t c + d, respectively. Here we use the max- 
imal positive gluon polarization of GRSV with Ag — g 
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Fig. 5. As in FigureObut now for the improved angular accep- 
tance of Smax ~ 180 mrad. 



at the input scale, for which the cancellation between 
the photon-gluon fusion and QCD Compton subprocesses 
mentioned above is less relevant. For our "default" gluon 
polarization used in Fig. (21 (b) the polarized cross sec- 
tion has a node at some pt,c such that ratios are difficult 
to visualize. From Fig. 01 one infers that the "direct" part 
dominates in absolute value in the entire pT,c-ra'nge shown 
and that the "resolved" and "direct" contributions have 
opposite signs. Turning to the individual subprocesses in 
the lower panel of Fig. ^ we note that the QCD Comp- 
ton 7g-channel always gives a positive contribution to the 
cross sectiou"'^ whereas the sign of the photon-gluon fusion 
channel is anti-correlated with the sign of Ag and its rel- 
evance scales with the magnitude of the unknown gluon 
polarization. This leads to a partial cancellation between 
the two "direct" channels for positive gluon polarizations 
which explains the smallness of the polarized cross sec- 
tion for the GRSV "standard" gluon observed in Fig. (21 
(b). Of the "resolved" processes only the scattering of a 
quark with large momentum fraction in the photon off 
a gluon in the nucleon makes a significant contribution, 
other channels are negligible. We note that in the unpo- 
larized case (not shown here) all contributions to the cross 
section are positive with the "direct" part accounting for 
80 percent or more. 

We also wish to comment on the momentum fractions 
Xb predominantly probed in the nucleon in a measurement 



^ dAatot < for the large and positive gluon polarization 
used in Fig.0]such that dAa-^q/ dAatot < 0. 
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Fig. 6. Scale dependence of the unpolarized (a) and polarized (b) photoproduction cross section for two charged hadrons in LO for 
Hermes. The scales in Eq. Q are varied simultaneously in the range {Pt,c + PT.d)/'^ < 4:{Pt,c + PT.d) indicated by the 



shaded bands. The solid lines correspond to the default choice where /i 
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of hadron-pair photoproduction at COMPASS. For the un- 
polarized cross section this can be easily specified by look- 
ing at the distribution in Xb for a given bin of pt,c- For 
example, ior pT,c around 1 GeV we find (xb) — 0.12 ±0.05 
which is consistent with the cc-range estimated by COM- 
PASS ^I]. However, similar estimates for the polarized 
cross section and hence for Ag/ g are impossible without 
knowing Ag beforehand as both the polarized cross sec- 
tions and the helicity parton distributions are not positive 
definite. The relevance of contributions of opposite sign to 
the cross section strongly depends on the gluon polariza- 
tion and completely obscures the meaning of an averaged 
{xb) here. This issue can be only consistently resolved in 
a future global analysis of polarized parton densities. 

Next it is interesting to check whether the upgrade 
of the angular acceptance of the Compass experiment to 
^max — 180 mrad will enhance their sensitivity to Ag. 
In Fig. [51 we present our expectations for the double-spin 
asymmetry All as a function of the transverse momentum 
Pt,c of one of the hadrons. Except for ^max = 180 mrad 
all other settings and cuts are the same as in Fig. |3| The 
statistical precision for such a measurement is again esti- 
mated with the help of Eq. H10|l for an integrated luminos- 
ity of 1 fb^^. From the upper panel of Fig.Oone infers that 
the "resolved" contribution modifies the asymmetry even 
more significantly than for 0,nax = 70 mrad. Even worse, 
there is now also a strong dependence on the model used 
to describe the parton content of the circularly polarized 
photon as can be seen in the lower panel of Fig.El This can 
be readily understood by noticing that due to the larger 
angular coverage one now probes the partonic structure 
of the photon also at momentum fractions x-y where the 
details of the unknown non-perturbative input do matter. 
Only in the high-p^ c-region, where x-^ — > 1, the depen- 
dence on the model for Ap becomes small. Clearly, a vi- 
able strategy would be to analyze data with different cuts 
on ^inax or for bins in and to learn as much as possible 
about Ag first. Data up to 6'max ~ 180 mrad might then 
be used for studying the non-perturbative structure of cir- 



cularly polarized photons. Needless to mention again, that 
the validity of the pQCD framework for two-hadron pro- 
duction at Compass has to be confirmed prior to studies 
of All- 



3.2 Two-Hadron Production at HERMES 

At the Hermes experiment at DESY longitudinally po- 
larized electrons/positrons with a beam energy of -Eg — 
27.5 GeV were scattered off both, a polarized deuterium 
or a polarized hydrogen gas target. The available c.m.s. 
energy of about ^/S ~ 7.5 GeV is lower than at COMPASS 
which even further limits the range of accessible trans- 
verse momenta. On average the lepton beam polarization 
is Ve — 53%. For the polarization of the gas target we 
take Vd ~ Vp ~ 85%, and, contrary to a sohd-state tar- 
get, there is no dilution, i.e., = Td = 1- 

We concentrate on phenomenological studies for a po- 
larized deuterium target in line with the data sample with 
the highest statistics in the Hermes spin physics program 
which came to an end recently. The Hermes experiment 
has an angular acceptance of 40 mrad < 0iab < 220 mrad 
for hadrons. For all our numerical studies we demand 
a transverse momentum of at least 1 GeV for both de- 
tected hadrons i?c,d- We choose a maximal photon vir- 
tuality of Q^ax = 0.1 GeV^ in Eq. and restrict y 
to 0.2 < y < 0.9. The fractions of the parent parton's 
momenta carried by the produced hadrons are Zc.d > 0.1 
Again, all scales in Eq. Q are set equal to 
unless stated otherwise. 

Figure|n|shows the dependence of both the unpolarized 
and polarized LO photoproduction cross section, Eq. Q, 
on the unphysical factorization/renormalization scales var- 
ied in the range {p'^^^ +PT,d)/^ < A*^ < '^{Pt,c + Pr.d)- 
Not unexpectedly, due to the smaller c.m.s. energy of the 
Hermes experiment, the scale dependence is even larger 
than for COMPASS, cf. Fig. El All remarks about poten- 
tial problems with the applicability of perturbative meth- 
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Fig. 7. As in Fig.Obut know for Hermes kinematics (see text). 
Estimates of the statistical accuracy are based on an integrated 
luminosity of 200 pb~^. 



ods at fixed-target energies and the need for unpolarized 
"benchmark" cross sections also apply here. 

Next we consider the corresponding double-spin asym- 
metry in Fig. 13 As in Fig. |2|we study the relevance 
of the "resolved" photon contribution in the upper panel 
and the dependence on models for the non-perturbative 
partonic structure of circularly polarized photons in the 
lower panel of Fig. [3 Estimates of the statistical accu- 
racy, Eq. (|1U|) . are based on the integrated luminosity of 
£ = 200 pb^^ actually collected by Hermes and the pa- 
rameters as specified above. Here, the "resolved" photon 
processes cause a much less pronounced shift in the asym- 
metry than for COMPASS, see Fig.O Also, there is almost 
no difference between the results obtained with the two 
extreme models for the Ap densities. This is because for 
the same transverse momentum "Pt.c the Hermes exper- 
iment is closer to the end of phase-space than Compass, 
i.e., xt.c = '^Pt,c/"/S is closer to one. On average Her- 
mes probes larger momentum fractions both in the nu- 
cleon and in the photon which explains our results. Our 
estimate of the statistical accuracy in Fig. [3 are such that 
only the region 1 < px^c < 1-5 GeV is potentially use- 
ful to further constrain the gluon polarization Ag. New 
experimental results from Hermes for both All and the 
underlying unpolarized cross section will become available 
in the near future ^Hl superseding the results of an earlier 
pubhcation [^j. 
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Fig. 8. As in Fig. O but now for Hermes kinematics (see text). 



We finish this section with a detailed study of the 
different contributions to the polarized photoproduction 
cross section. The upper and lower panel of Fig.|Hlshow the 
fractional contributions of dAadiT and dAaics and of the 
different partonic LO channels a + b c + d, respectively. 
As for Fig. 01 we use the maximal positive gluon polariza- 
tion of GRSV with Ag = g a.t the input scale. Again the 
choice is due to the fact that for our "default" gluon po- 
larization, GRSV "standard" , the cross section develops 
a node at some pt,c such that ratios are difficult to visu- 
alize. As is already expected from Fig. |3 the "resolved" 
photon cross section is fairly small throughout, though 
not completely negligible. Again it contributes with the 
opposite sign than the "direct" cross section, cf. Fig. 0] 
As for Compass kinematics before, we observe a partial 
cancellation of the two "direct" channels: photon-gluon 
fusion and QCD Compton scattering. None of the four 
"resolved" contributions makes a significant contribution 
at LO. 



3.3 Prospects for a Future Polarized ep Collider 

The most interesting option for a future experimental spin 
physics facility is a first polarized lepton-proton collider 
such as the eRHIC project at BNL 15 currently under 
discussion. Here we consider the asymmetric collider op- 
tion using the existing 250 GeV proton beam of RHIC 
and a new 10 GeV electron beam, i.e., a c.m.s. energy of 
y/S = 100 GeV. The physics program of such a machine is 
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Fig. 9. Double-spin asymmetry for photoproduction of two 
charged hadrons in LO at \fS — 100 GeV. Results are given for 
two different gluon polarizations of the nucleon and the two sce- 
narios for the Af densities. Both hadrons are required to have a 
transverse momentum of at least 2 GeV. The laboratory rapidity 
of one of the hadrons is integrated in the range 0.6 < rjd < 2.6 
and 2.6 < rjd < 3.6 in the upper and lower panel, respectively. 
Estimates of the statistical accuracy refer to an integrated lumi- 
nosity of 1 fb^^. 



ger" hadron is integrated in the range 0.6 < rjd < 2.6 and 
2-6 < rjii < 3.6 in the upper and lower panel, respectively, 
and the rapidity of the other hadron is left unintegrated. 
For ?7c ^ our results are fairly independent of the choice 
of a particular model for Af while for r/c ^ 1 the results 
strongly depend on Z\/^. This is in particular the case if 
the trigger hadron is detected more forward in the incom- 
ing proton's direction, i.e., for 2.6 < rjd < 3.6 (lower panel 
of Fig. ISJ . Note that in this subsection positive rapidities 
denote the proton direction, in line with common conven- 
tions usually used at HERA. Again this behavior can be 
understood by looking at the typical momentum fractions 

probed in Eq. (jSJ: for rjc^d large and positive one finds 

^ 1 and for rjc^d large and negative approaches one. 

We also give estimates for the statistical accuracy in 
Fig. 121 based on beam polarizations of Ve,p — 0.7 and 
an integrated luminosity of 1 fb~^. The latter is expected 
to be accumulated within only a few weeks of running 
eRHIC so that the statistical accuracy can be eventually 
much better than the one shown in Fig. |51 In addition, 
we demand a minimum transverse momentum for both 
hadrons of 2 GeV and Zc.d > 0.1. For the equivalent pho- 
ton spectrum in Eq. ^ we use similar parameters as the 
HI and ZEUS experiments at HERA: Q^^^ = 0.5 GeV^ 
and the momentum fraction taken by the photon is limited 
to 0.2 < y < 0.85. Note that the "boost" factor between 
the laboratory and the ep-c.m.s. frame is very similar for 
the asymmetric collider option for eRHIC and for HERA. 

Keeping in mind that the polarized gluon distribution 
Ag should be known fairly well in the range relevant in 
Fig.El i.e. X > 0.01, from RHIC by the time eRHIC would 
start to operate, there are excellent prospects to study 
the so far unknown parton content of circularly polarized 
photons. 



based on the extensive and highly successful exploration of 
unpolarized ep collisions at the DESY-HERA collider. The 
determination of the gluon density from scaling violations 
in DIS down to very small x and establishing the concept 
of photonic parton densities in photoproduction processes 
are some of the many physics highlights of HERA. 

From similar precision studies of scaling violations in 
polarized DIS at eRHIC it would be possible to access the 
gluon polarization down to x ~ 10^"^, about one decade in 
a; lower than in pp collisions at RHIC. In photoproduction 
processes it is expected that "resolved" contributions play 
a much more significant role at eRHIC than at fixed-target 
experiments jl6| . This offers exciting prospects to learn 
about the Af^ densities as was already demonstrated in 
case of single-inclusive pion photoproduction in Ref. 

Studies of hadron-pair production have the advantage 
that one has a better control on the momentum fractions 
probed in the hadron and in the photon. For instance, by 
demanding a "trigger hadron" in the proton's forward di- 
rection and scanning the other hadron's rapidity one can 
select kinematical regions which are particularly sensitive 
to the non-perturbative structure of the photon. This is 
demonstrated in Fig. |51 where we give estimates for All 
at VS* = 100 GeV. The laboratory rapidity of the "trig- 



4 Conclusions 

In summary, we have presented a phenomenological study 
of spin-dependent photoproduction of hadron pairs at 
c.m.s. energies relevant for the Compass and Hermes 
experiments as well as for a possible future polarized ep 
collider. So far, our studies are limited to the LO approx- 
imation of QCD but will be amended to NLO accuracy in 
due time. 

We have consistently included the "direct" and "re- 
solved" photon contributions to the photoproduction cross 
section. It turned out that the "resolved" part leads to 
a significant shift in the experimentally relevant double- 
spin asymmetries and has to be accounted for in future 
analyzes of data. Fixed-target experiments are, however, 
mainly sensitive to the perturbative "pointlike" part of the 
photon structure. This simplifies attempts to extract the 
gluon polarization Ag which is the main goal of COMPASS 
and one of the goals of Hermes. The non-perturbative 
parton content of circularly polarized photons can be pro- 
bed in detail in photoproduction processes at higher c.m.s. 
energies which hopefully will become available at some 
point in the future. By that time we most likely have a 
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good knowledge of the spin structure of the nucleon, in 
particular, from ongoing measurements at RHIC. 

The double-spin asymmetries for both Compass and 
Hermes show the expected sensitivity to the gluon polar- 
ization in the nucleon. Keeping in mind the very sizable 
scale dependence at LO all results for All have to be 
taken with a grain of salt unless the applicability of per- 
turbative methods for two-hadron production at compar- 
atively low transverse momenta and c.m.s. energies has 
been thoroughly investigated and demonstrated. This is 
best achieved by comparing the underlying unpolarized 
cross sections for the production of hadron-pairs with the- 
oretical expectations. If these checks are passed, the pro- 
duction of hadron pairs in lepton-nucleon collisions will 
be an interesting and complementary tool to further con- 
strain the polarized gluon density at momentum fractions 
of about X = 0.1-^0.2. If it turns out that data and theory 
do not match, these measurements will open up a window 
to study the effects of all-order resummations, the rele- 
vance of higher- twist corrections, and perhaps the transi- 
tion to the non-perturbative regime so far little explored 
and understood. 
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